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On the atmospheric fragmentation of small asteroids 
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Abstract. It is known, from observational data recorded from airbursts, that small asteroids breakup at dynamical 
pressures lower than their mechanical strength. This means that actual theoretical models are inconsistent with 
observations. In this paper, we present a detailed discussion about data recorded from airbursts and about several 
theoretical models. We extend and improve a theory previously outlined for the fragmentation of small asteroids 
in the Earth atmosphere. The new condition for fragmentation is given by the shock wave-turbulence interaction, 
which results in sudden outburst of the dynamical pressure. 
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1. Introduction 

Collisions between cosmic bodies are divided into two 
regimes: gravity- and strength-dominated regime. But the 
transition from one regime to the other is not well known 
(see, for example, Fig. 5 in Durda et al. 1998). In the case 
of collision with Earth, the presence of the atmosphere 
makes things more difficult. Observations of airburst of 
small asteroids (up to tens of metres) show that the frag- 
mentation occurs when the dynamical pressure is lower 
than the mechanical strength and this conundrum has 
not a sa tisfac tory explanation yet (see Ceplecha 1996b , 
Foschini 2000). It is worth noting that also the fracture 
itself has still many unknown features (for a review, see 
Fineberg & Marder 1999). 

Studies on the fragmentation of small asteroids have 
also a great importance in impact hazard. Although the 
damage caused by Tunguska-like objects can be defined as 
"local" , it is not negligible. The Tunguska event of 30 June 
1908 resulted in the devastation of an area of 2150±25 km 2 
and the destruction of more than 80 million trees (for a 
review, see Trayner 1997, Vasilyev 1998). Still today there 
is a wide debate all over the world about the nature of 
the cosmic body which caused that disaster. Just on July 
1999 an Italian scientific expedition, Tunguska99, went to 
Siberia to collect data and samples (Longo et al. 1999 ). 

In this article, we discuss firstly the current models 
on atmospheric fragmentation of small asteroids and in- 
consistencies with observations (Sect. 2). In Sect. 3, we 
deal with the question of strength and aerodynamic load. 
In Sect. 4, we extend the approach previously outlined 
in Foschini ( |1999h| ), thereafter called Paper I. We study 



the hypersonic flow around a small asteroid, with particu- 
lar reference to the definition of the pressure crushing on a 
cosmic body (Sect. 5). Therefore (Sect. 6), it is possible to 
obtain the condition for fragmentation under steady state 
conditions, that was already outlined in the Paper I. The 
Sect. 7 deals with the turbulence and massive ablation, 
putting the basis for the condition of fragmentation under 
unsteady conditions (Sect. 8). Some numerical examples 
conclude the paper (Sect. 9). 

2. Problems with current models 

Present models consider that the fragmentation begins 
when the condition: 



PooV- = S 



(1) 
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is satisfied. In Eq. (H), is the density of undisturbed 
stream, V is the speed of the body, S is the material me- 
chanical strength. The term PooV 2 is referred as the dy- 
namical pressure in the front of the cosmic body. 

However, observations of very bright bolides prove 
that large meteoroids or small asteroids breakup at dy- 
namical pressures lower than their mechanical strength. 
Today there is still no explanation for this conundrum. 
This is a scientific problem of great interest, but it is also 
of paramount importance, because it allows us to know 
whether or not an asteroid might reach the Earth surface. 
In addition to this, the atmospheric breakup has also ef- 
fect on the crater field formation (Passey & Melosh 1980) 
or on the area devastated by the air blast. 

The interaction of a cosmic body in the Earth atmo- 
sphere in the strength-dominated regime can be divided 
into two parts, according to the body dimensions. For 
millimetre- to metre-sized bodies (meteoroids), the most 
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Table 1. Meteoroid strength category. After Ceplecha et 
al. (1993) 



Category 


Range of p x V z [MPa] 


Mean p^V [MPa] 


a 


Pao V' 2 < 0.14 


0.08 


b 


0.14 < Poo V 2 < 0.39 


0.25 


c 


0.39 < p^V 2 < 0.67 


0.53 


el 


0.67 < p^V 2 < 0.97 


0.80 


e 


0.97 < pooV 2 < 1.20 


1.10 



useful theoretical model is the gross-fragmentation model 
developed by Ceplecha et al. (1993) and Ceplecha (1999). 
In this model, there are two basic fragmentation phenom- 
ena: continuous fragmentation, which is the main process 
of the meteoroid ablation, and sudden fragmentation or 
the discrete fragmentation at a certain point. 

For small asteroids another model is used, where the 
ablation is contained in the form of explosive fragmenta- 
tion, while at high atmospheric heights it is considered 
negligible. Several models have been developed: Baldwin 
& Shaeffer (1971), Grigoryan (1979), Chyba et al. (1993), 
Hills & Goda (1993), Lyne et al. (1996). A comparative 
study on models by Grigoryan, Hills & Goda, and Chyba- 
Thomas-Zahnle was carried out by Bronshten (1995). He 
notes that the model proposed by Chyba et al. does 
not take into account fragmentation: therefore, the de- 
struction heights are overestimated (about 10-12 km). 
Bronshten also concludes that Grigoryan and Hills-Goda's 
models are equivalent. 

Despite the particular features of each model, fragmen- 
tation is always considered to start when Eq. ([I]) is satis- 
fied. 

Although direct observations of asteroid impact are 
not available, it is possible to compare these models with 
observations of bodies with dimensions of several me- 
tres or tens of metres. Indeed, in this range, the gross- 
fragmentation model overlaps the explosive fragmentation 
models. As underlined several times by Ceplecha (1994, 
1995, 1996b), observations clearly show that meteoroids 
breakup at dynamical pressures lower (10 times and more) 
than their mechanical strength. These data are obtained 
from photographic observation of bright bolides and the 
application of the gross-fragmentation model, that can be 
very precise. 

According to Ceplecha et al. (1993) it is possible to 
distinguish five strength categories with an average dy- 
namical pressure of fragmentation (Table [l]) . For contin- 
uous fragmentation the results obtained also indicate that 
the maximum dynamical pressure is below 1.2 MPa. Five 
exceptions were found: 4 bolides reached 1.5 MPa and one 
survived up to 5 MPa (Ceplecha et al. 1993). 

It would be also very important to relate the ablation 
coefficient with the dynamical pressure pocV 2 at the frag- 
mentation point, in order to find a relationship between 
the meteoroid composition and its resistance to the air 
flow. To our knowledge, a detailed statistical analysis on 
this subject does not exist, but in the paper by Ceplecha 



et al. (1993) we can find a plot made by considering data 
on 30 bolides (we refer to Fig. 12 in that paper). We note 
that stony bodies (type I) have a wide range of PooV 2 val- 
ues at the fragmentation point. In the case of weak bodies, 
we can see that there is only one cometary bolide (type 
IIIA), but this is due to two factors: firstly, cometary bod- 
ies undergo continuous fragmentation, rather than a dis- 
crete breakup at certain points. Therefore, it is incorrect 
to speak about fragmentation pressure; we should use the 
maximum tolerable pressure. The second reason is that 
there is a selection effect. Indeed, from statistical studies, 
Ceplecha et al. (1997) found that a large part of bod- 
ies in the size range from 2 to 15 m are weak cometary 
bodies. However, a recent paper has shown that statistics 
from physical properties can lead to different results when 
compared with statistics from orbital evolution (Foschini 
et al. 2000). To be more precise, physical parameters prove 
that, as indicated above, a large part of small near Earth 
objects are weak cometary bodies, whilst the analysis of 
orbital evolution proves a strong asteroidal component. 

In addition to data published in the paper by Ceplecha 
et al. (1993) and Ceplecha (1994) we can see in Table || 
some specific cases of bright bolide s; for details, see the 
papers quoted or to Foschini ( [2000 ). 



3. Stresses and strengths 

For the sake of the simplicity, in the above section, we 
adopted commonly used values for mechanical strength, 
i.e. 1 MPa for cometary bodies, 10 MPa for carbonaceous 
chondrites, 50 MPa for stony bodies, and 200 MPa for 



iron bodies (see, for example, Hills & Goda 1993 ). Only 
Bronshten (200C) rejected the values for cometary bodies, 
proposing the ra nge 0.0 2 — 0.4 MPa. He used values cal- 
culated by Opik ( 1966 ) after observation of tidal disrup- 
tion of Sun-grazing comets under the gravity of the Sun. 
However, it is worth noting that the strength obeys to scal- 
ing laws: therefore the larger is the body, the smaller is 
the strength. Opik's calculations refer to the comet Ikeya- 
Scki (1965f), which has an estimated nucleus, before the 
breakup, of about 8.3 km. Small asteroids are in the range 
of several tens of metres, up to some hundreds of metres. 
The Tunguska Cosmic Body (TCB), to which Bronsthcn 
apply his calculations, is in the range 50 — 100 m. We 
can calculate how the scaling law changes the value of the 
strength. The formula is shown in the paper by Tsvetkov 



Skripnik (1991) and we recall it for simplicity: 



S = S' 



(2) 



where m and m! are the masses of the bodies and a is 
the scale factor: the more inho mogen eous is the mate- 
rial, the higher is a. Turcotte ( 1986 ) adopted a value 
of a = 0.12 for glacial tills and we use this value. The 
comet mass can be calculated as 3 • 10 14 kg, by using a 
density of 1000 kg/m 3 , while the Tunguska Cosmic Bod y 
has an estimated value around 10 s kg (Vasilyev 1998 ). 
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Name 


Date 


Tr? 

PooV 


s 


Source 


Pribram 


Apr 7, 1959 


9.2 


50 


ReVelle (1979) 


Lost City 


Jan 3, 1970 


1.5 


50 


Ceplecha (1996a) 


Sumava 


Dec 4, 1974 


0.14 


1 


Borovicka & Spurny (1996) 


Innisfree 


Feb 6, 1977 


1.8 


10 


ReVelle (1979) 


Space based obs. 


Apr 15, 1988 


2.0 


50 


Nemtchinov et al. (1997) 


Space based obs. 


Oct 1, 1990 


1.5 


50 


Nemtchinov et al. (1997) 


Benesov 


May 7, 1991 


0.5 


10 


Borovicka et al. (1998a, b) 


Peekskill 


Oct 9, 1992 


1.0 


30 


Ceplecha et al. (1996) 


Marshall Isl. 


Feb 1, 1994 


15 


200 


Nemtchinov et al. (1997) 



With these values in Eq. (|2|) we obtain a value of the me- 
chanical strength for a Tunguska-sized body, in the range 
0.12 — 2.4 MPa, that is compatible with commonly used 
value of 1 MPa. 

The scaling law for mechanical strength derives from 
the assumption that the fragmentation was a process of 
consecutive elimination of defects under increasing load. 
Baldwin & Sheaffer (1971) consider that the reason for 
the presence of cosmic bodies with very low fragmenta- 
tion pressure can be explained by the assumption that 
additional flaws and cracks may be created by collisions 
in space, even though they do not completely destroy the 
cosmic body. 



Tsvetkov & Skripnik (1991) made an interesting study 
on the fragmentation according to strength theory and 
scaling laws, but this kind of study is useful only for mete- 
orites and cannot explain airbursts. They also show that 
the condition of fragmentation S — PooV 2 is not valid 
even under the assumptions of scaling laws. Indeed, they 
showed that the aerodynamic loading never reach the ul- 
timate strength. For this reason, they searched the cause 
of fragmentation in the particular structure and extreme 
inhomogeneity of cosmic bodies. With an appropriate se- 
lection of the scale factor, Tsvtkov & Skripnik obtained a 
good fit. 

However, the shock compression and heating during 
the atmospheric entry will result in an elimination of inter- 
nal cracks, making the body more compact (see, for exam- 
ple, Davison & Graham 1979 or Zel'dovich & Raizer 1966). 
If internal voids are so large to survive to shock compres- 
sion, they could give an explanation for some episodes, 
but not a general theory (Foschini 1998). 

4. The hypersonic flow approach 

Almost all models described in Sect. 2 deal with the mo- 
tion of a cosmic body in the Earth atmosphere. However, 
it is worth noting that we cannot observe directly the cos- 
mic body: we can only see the light emitted during its 
atmospheric entry. Therefore, we have to introduce in the 
equations several coefficients, that cannot be derived di- 
rectly from observations. 

If we turn our attention to the hypersonic flow around 
the body, we could have data from direct observations. 
Among models discussed above, only Nemtchinov et al. 



(1997) tried to investigate the hypersonic flow around the 
asteroid with a numerical model. Foschini (1999b) inves- 
tigated the analytic approach: indeed, although the de- 
tails of hypersonic flow are very difficult to calculate and 
there is need of numerical models, the pressure distribu- 
tion can be evaluated with reasonable precision by means 
of approximate methods. In the limit of a strong shock 
(M >> 1) several equations tend to asymptotic values 
and calculations become easier. For example, the ratio of 
densities across the shock is: 



Po 

Poo 



7 + 1 
1 



7 



(3) 



where po is the density in the stagnation point and 7 is 
the specific heats ratio. 

In the Paper I, it was showed the crucial role of the 
temperature, instead of the pressure, in the stagnation 
point. The gas is in local thermodynamic equilibrium 
(LTE), that is, matter is in equilibrium with itself, but 
not with radiation, which can escape. Particle densities 
depend on the temperature only and it is possible to use 
Boltzmann or Saha equations. The temperature is calcu- 
lated directly from observations of spectra (e.g. Borovicka 
& Spurny 1996, Borovicka et al. 1998b) and, thus, it is a 
reliable starting point. 

Particularly, as stated in the Paper I, the temperature 
at the stagnation point is very important. Changes in the 
stream properties are mainly due to changes in the stagna- 
tion temperature, which is a direct measure of the amount 
of the heat transfer. The enthalpy change Ah is: 



Ah = c P AT 



(4) 



where cp is the specific heat with constant pressure. From 
Eq. (||) , it is possible to relate the maximum speed of the 
stream (which is close to the cosmic body speed) to the 
stagnation temperature (see Paper I): 



V m 



27 
7-1 



RT 



(5) 



where R is the constant of the specific gas. The Eq. (||) is 
valid for an adiabatic flow whether reversible or not. 

Under certain approximations, the flow can be con- 
sidered adiabatic: heat transfer in hypersonic motion is 
strongly reduced. The boundary layer, where convective 
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heat transfer takes place, is very thin and can be consid- 
ered negligible. 

On the other hand, if we consider the Eq. (||), we see 
that it is satisfied for unrealistic high temperatures. For 
example, for a speed of V = 12 km/s in air we obtain 
To = 72000 K. The introduction of ablated matter make 
7 higher and R lower, resulting in a further increase of To. 

However, at high temperatures some effects of the gas 
become important: the main effects are the dissociation 
and ionization of air molecules. For 2000 < T < 4000 K 
the oxygen molecules break down to single atoms; for 
4000 < T < 9000 K there is the dissociation of N 2 ; 
for T > 9000 K ionization occurs (see, for example, 
Oosthuizen & Carscallen 1997). In the Ah term of the 
Eq. (||), we should take into account also the heat ab- 
sorbed by these processes. As a result, To drops to more 
reasonable 10000 - 15000 K or so. 

The treatment of real gas properties is probably one 
of the most difficult problem in aerothermodynamics and 
involves complex numerical modelling (for a review, see 
Tirsky 1993, Gnoffo 1999). 

5. The equation of state 

Before going on, it is necessary to evaluate the state of 
the shocked gas, in order to select the more appropriate 
equation of state. 

From the point of view of pressure, we can see that 
the gas around an asteroid during the atmospheric entry, is 
still an ideal gas. Indeed, the limit of pressure to become a 
degenerate gas is around 10 14 Pa (Eliezer 1991). A simple 
calculation shows that the asteroid reaches the maximum 
dynamical pressure at the sea level (p s \ = 1.293 kg/m 3 ) 
and if it has the maximum geocentric speed (V = 72 km/s, 
although it is quite impossible to have an asteroid with 
such a speed). Therefore, the maximum PooV 2 « 7T0 9 Pa, 
about four orders of magnitude below the limit of degen- 
eracy. 

From the point of view of the temperature, the ques- 
tion is a bit difficult to handle. From spectroscopic ob- 
servations of bright bolides (Borovicka & Spurny 1996, 
Borovicka et al. 1998b), we can see that the spectrum is 
composed by two temperatures emission: the first value 
lies in the range between 4000 — 6000 K and the second 
value is about 10000 K. The temperature of 10000 K is 
about the limit between normal and partially ionized gas. 
Therefore, we shall adopt a more general point of view, 
based on microscopic physics. 

From the kinetic theory, we know that the pressure is 
additive, therefore we can separate it in contributions from 
different species composing the gas: electrons, photons, 
ions, and atoms. We can write: 



P = Pi + P c + Pr 



j 



AjfcT +N c kT c + 



aT 4 



(6) 



where k is the Boltzmann constant and a is the radia- 
tion constant, related to the Stefan Boltzmann constant 



(a = Aa/c = 7.56591 • 10~ 16 [Pa • K~ 4 ]); N is the volume 
density of species (the index e is for electrons); the sum 
over the index j concerns all heavy particles (ions and 
neutral atoms). 

Eq. (§) has been obtained by means of quite general 
assumptions, without any particular restrictions (more de- 
tails can be found in many places; for example, in Lang 
1999). We can apply Eq. (j^) to hypersonic flow around an 
asteroid in the Earth atmosphere. The only restrictions are 
determined by the momentum distributions of particles: 
ions, atoms, and electrons obey generally to a Maxwellian 
distribution, photons obey to the Bose-Einstein statistics. 
These conditions are satisfied in processes occurring dur- 
ing the atmospheric entry of a cosmic body. 

6. The pressure of fragmentation 

We can note that, for typical temperatures in airbursts, 
the radiation pressure is negligible. For T = 10000 K, 
we have that P r « 2.5 Pa, which is negligible when com- 
pared to PooV 2 at the fragmentation point (see Sect. 2). 
Therefore, for airbursts and meteors, the Eq. (^) becomes: 



P » ^ + P c 



(7) 



In order to calculate the densities of species, we have 
to take into account that the fluid around the asteroid is 
in local thermodynamic equilibrium (LTE) and we can use 
the Saha relation: 



N c Ni _ (27rm£T) 3 / 2 2 gi eE { 



(8) 



where 2V a is the density of neutral atoms, E\ is the ion- 
ization potential, g a and gi are respectively the statistical 
weight of the ground state of the neutral atom and of the 
ion. Some of these values are listed in Table ||. 

Table 3. Values of E- u g a , and for some species in the 
flow 



Species 


Ei [eV] 




m 


Na 


5.14 


2 


l 


K 


4.34 


2 


i 


0, 


12.05 


3 


4 


N 2 


15.6 


1 


2 



For each airburst, with a given composition, we can 
calculate densities of all species and obtain the effective 
pressure. But this is not the purpose of this work. We will 
introduce some hypoteses for the sake of the simplicity, in 
order to do some calculations of order of magnitude. 

The Saha equation (^|) can be rewritten in terms of 
degree of ionization a: 



N tot a 2 = (2^mfcT) 3 / 2 2 9i 
(1 - a) W g a 



exp(- 



eEi 
~kT 



(9) 




where N tot is the total number density of particles; there- 
fore, N c — aN tot and so on. The Eq. (]?]) can be rewritten: 



P = (1 + a)N tot kT = (1 + a)pRT 



(10) 



that is the well known equation of state for a gas with a 
degree of ionization a. 

Foschini (1999a) underlined the importance of alkaline 
metals in the production of electrons in meteors. In addi- 
tion, von Zahn et al. ( |1999 ) recently showed the particular 
role of the potassium in Leonid meteors. Indeed, alkaline 
metals, owing to their low ionization potential easily ionize 
even at low temperatures. For example, we assume that 
JVtot = 6.022 • 10 23 m~ 3 (it is the Avogadro number) and 
we calculate the degree of ionization of a gas composed 
by single species listed in Table g. For a temperature 
T = 10000 K, we have: for sodium, a = 0.92; for potas- 
sium, a — 0.96; for oxygen, a = 0.1 and for nitrogen, 
a = 0.01. 

It is known that, in meteors, the impact ionization 
from impinging molecules is very important. A nitrogen 
molecule with a speed of 18 km/s has an energy of about 
45 eV. However, as the asteroid penetrates deeper in the 
atmosphere, the flow generates a shield that protects the 
asteroid from direct impacts (for a detailed description see 
Ceplecha et al. 1998). Therefore, the thermal ionization 
become important particularly in small asteroids/comets, 
when impact ionization is strongly reduced. So, alkaline 
metals play an important role, because they provide a 
huge amount of electrons. 



We have also to take into account that we are dealing 
with particular condition in the stagnation point, where 
there is the maximum thermo-mechanical stress. In this 
point, we can consider the fluid as highly ionized: the ab- 
lation of alkaline metals provide a source of ions for the 
flow. In the stagnation point - not in the whole flow - we 
can consider the gas with a de gree of ionization a — ► 1. 
Therefore, by substituting Eq. ( jlp| ) in Eq. ([|), we obtain 
the condition for fragmentation: 



Mtiax — 



20 



P 



7 — 1 (1 + a)p 



(11) 



Eq. ([y]) was already outlined in the Paper I, where it 
was considered a = 1. 

It is worth noting that Eq. (11) is the condition of 



fragmentation, not the condition for the airburst, as de- 
duced by Bronshtcn ( 200C ). The airburst generally occurs 
after a scale height, as shown by several studies on super- 
bolides. In the Paper I, there was an unfortunate error in 
considering the explosion height, instead of the fragmen- 
tation height (although it has negligible effect, as we shall 
see in the Sect. 9). Therefore, Eq. flll| ) must be revised as 
follows: 



Knax — 



20 



7-1(1 + a)p s \ 



exp( 



h + H 
H ] 



(12) 



where p s \ is the air density at sea level, h is the airburst 
height, and H is the scale height. 
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Before to analyse the Eq. (|l2|), it is necessary to look 
at the value for 7. In the Paper I, 7 = 1.7 was used and 
it was taken from experimental measurement in hyperve- 
locity impacts. Bronshten (2000) does not agree with this 
value and invokes argument to support a value of 7 = 1.15, 
that is compatible with air at high temperatures and un- 
der shock loading. Bronshten's assumption are reasonable 
when we have a diatomic gas at high temperature: the 
rotation and vibration of air molecules change the val- 
ues for 7. But, if we have a monatomic gas or metal va- 
pors, we have to consider a 7 = 5/3, unless the temper- 
ature is so high to change ap preci ably the electron en- 
ergy (see Zel'dovich & Raizier |l966 ); however, this is not 
our case. This is confirmed also by meteor spectra: they 
show that the surrounding gas is monatomic, although 
mole cules bands are sometime recorded (Ceplecha et al. 



199S 



It is worth noting that the question about the value 
of 7 is not simple. Perhaps, it will be solved when experi- 
mental data will be available. 

7. Turbulence and massive ablation 

Eq. (|lj) shows that it is necessary to have higher speeds 
in order to reach the value of the mechanical strength. It 
is worth noting that Eq. (12) is valid under steady state 
condition. In this case, we can neglect the contribution 
of turbulence, because of compressibility effects at high 
Mach number (see Andreopulos et al. 20001 ). 

The influence of the turbulence in the large meteoroid 
entry was first analyzed by ReVelle ( 1979] ). In his ap- 
proach based on the motion of a single body in the at- 
mosphere, ReVelle studied how the convective heat trans- 
fer depends on turbulence. He found that the turbulent 
convective layer is negligible, except for speed lower than 
about 20 km/s. 

However, the pr esence of massive ablation changes the 
flow. Gupta ( 1983 ) noted that there are few experimen- 
tal data available and, most important, they are obtained 
with freestream Mach number in the range 3-7 and with 
negligible ablation. There are less data about experiments 
with moderate ablation, but with lower Mach number (up 
to 2.6). Starting from these data, models predict that con- 
vective heating is negligible with almost any turbulence 
model. This was the situation in 1983, but more recent re- 



views (Gnoffo 1999) show that there are still several things 
unknown. The only news is the entry of the Galileo probe 
in the Jovian atmosphere, which gave us useful informa- 
tion about hypersonic motion with massive ablation. The 
probe entered in the atmosphere at a relative speed of 
47.5 km/s and experienced an ablation rate of 7.4 km/s, 
with a total mass loss of about 79 kg (Gnoffo |l999| ). 

The probe did not suffer any fragmentation, despite of 
massive ablation, although it is still not clear the role of all 
processes during the entry. The layer between the shock 
and the probe is expected to be turbulent over almost the 
entire length of the Galileo probe, owing to the massiv e 
ablation and the large Reynolds number (Gupta 1983|) . 



The interaction of shock waves with turbulence leads to 
amplification of speed fluctuations and changes in length 
scales. It would be better to say that this occur when the 
flux is unsteady and the shock wave is subject to strong 
distortions. At high Mach number, but steady state mo- 
tion, the effect of compressibility takes place and there is 
no amplification; we can adopt Eq. (fP2]). 

But the atmospheric motion is generally unsteady. 
Fig. [I] shows the air temperature from tables of the US 
Standard Atmosphere 1976. The air temperature is instru- 
mental to derive the local sound speed (a — s/jBT), as a 
function of the height, for different latitude and seasons. 
As the temperature changes, the sound speed changes. 
The US Standard Atmosphere gives also a reference at- 
mosphere - and therefore a reference temperature - from 
which we can calculate a reference sound speed (Fig. ||). 

It is known that, for a small asteroid, the decelera- 
tion due to ablation is negligible; therefore we can con- 
sider that the body speed does not change until the frag- 
mentation begins. For a given asteroid, the Mach num- 
ber depends only on the local sound speed. In Fig. |^ 
we can see the Mach number for different values of the 
speed of the cosmic body. It results that the Mach number 
change substantially for higher asteroid speed, while for 
lower speed, changes of M become smooth. We note also 
strong variations when the body crosses the mesopause, 
the stratopause, and the tropopause. 



Ceplecha (1994) gave the mean values of end heights 



for large meteoroids (up to about 7 m): he found 32 km for 
type I; 43 km for type II; 57 km for type IIIA; and 69 km 
for type IIIB. These are end heights and if we consider that 
the fragmentation begins about a scale height above, we 
infer that the mean fragmentation height is: 39 km for type 
I; 50 km for type II; 64 km for type IIIA; and 77 km for 
type IIIB. That is, we can divide these bodies into two cat- 
egories: a first category, made with type II and type IIIB, 
which contains bodies that breakup during the crossing 
of the stratopause and the mesopause, respectively. The 
second category, made with type I and type IIIA, which 
contains bodies that breakup with a certain delay after 
the crossing of atmospheric pauses. The delay can be ex- 
plained with the different ability to ablate, and therefore 
to produce the turbulent boundary layer. Indeed, for bod- 
ies that breakup at the mesosphere (cometary bodies), the 
type IIIB has an average ablation coefficient higher than 
the type IIIA. For asteroidal bodies, that breakup at the 
stratosphere, the type II has a higher ablation coefficient 
than the type I. We should also take into account that the 
transition from laminar to turbulent flow depends on the 
body dimensions and speed. Large bodies at high speed 



reach the transition at higher heights (ReVelle 1979). 

It is worth noting that there are also other factors that 
can change the fragmentation height, such as the rotation 
(Adolfsson & Gustafson |1994| , Ceplecha |l996aD , but they 
are independent from the body type and can explain only 
deviations from the average values. 
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Fig. 2. Local sound speed from the reference temperature of the US Standard Atmosphere 1976. 
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Fig. 3. Mach number as a function of the height for different asteroid speed. Ml=12 km/s; M2=40 km/s; M3—72 km/s. 
8. The condition for fragmentation 



To find a quantitative condition for fragmentation by tak- 
ing into account the turbulence, at least for first order cal- 



culations, is very difficult. Problems derive from the lack of 
experimental data about the shock wave-turbulence inter- 



8 



Luigi Foschini: On the atmospheric fragmentation of small asteroids 



action at very high Mach number. Moreover, the few fea- 
sible experiments at moderate Mach number are strongly 
dependent on measurement systems and available data 
are insufficient to fully describe the mutual interaction 
of shock waves and turbulence. Last, but not least, the 
turbulence itself is one of the oldest unsolved problems in 
the history of science. In recent years, numerical models 
allowed detailed investigation, but the problem is com- 
plicated by the fact that, when dealing with turbulence, 
the averaging of governing equations introduces new un- 
knowns. Therefore, the number of available equations is 
not sufficient and it is necessary to assume a closure con- 
dition. The lacking of experimental data make it hard to 
make hypoteses on the closure condition and therefore 
numerical models are often contradict ory. F or revi ews o f 
these problems see An drcop ulos et al. ( 2000|) , Lele ( |l994|) , 
Adamson & Messiter ( 1980 ), and references therein. 

Despite of differences, it is clear that the shock wave- 
turbulence interaction produce an amplification of fluctu- 
ations. The amplification depends on the shock strength, 
the state of the turbulence, and its level of compressibil- 
ity. The most important outcomes are the amplification 
of velocity fluctuations and changes in the length scales 
(Andreopulos et al. 2000| ). This leads to changes in the 
dynamical pressure in the front of the asteroid, but also 
changes in pressure along the flank of the cosmic body; 
these changes can be further amplified by local irregular- 
ity, such as small "hills" . 

However, for the sake of simplicity, we make the as- 
sumption that the zone of greatest stress is the stagnation 
point. Or, it would be better to say, that we will continue 
to hold this hypotesis. Being under unsteady conditions we 
cannot apply isentropic relations, i.e. Eq. (ft2l). We start 
from the value of the pressure in the shock layer derived 
from the Rankine-Hugoniot relations in the hypersonic 
limit (P PooV 2 ). Being the gas cap around the asteroid 
a plasma, we have to consider the Eq. (|Io|), so we have: 



(13) 



For the sake of the simplicity, we assume that turbu- 
lence does not affect density; so we can consider that shock 
wave-turbulence interaction does affect only velocity. The 
problem is to set up the value of this amplification: we 
then introduce an amplification factor k to evaluate: 



P « (1 + a) Poc KV 2 



(14) 



As written above, we have no experimental data and 
numerical models are often contradictory. We can try to 
set up lower and upper limits. Rotman ( |1991[ ) reports am- 
plification of the kinetic energy of about 2 — 2.15. Jacquin 
et al. ( |1993| ) found that the amplification of the kinetic en- 
ergy depends on the density ratio; the factor can be up to 
12.7 for diatomic gases. For monatomic gases and plasmas 
the upper limit of the amplification is 6. 

In conclusion, we can assume for kinetic energy 2 < 
k < 6 and that this amplification value is valid also for 
pressure. Therefore, under distortion of the shock wave, 



the turbulence can leads to amplification of dynamical 
pressure up to 12 times the nominal value for a neutral 
gas (we have taken into account also the multiplicative 
factor 1 + a). Comparing with experimental data of the 
fragmentation of asteroids showed in Sect. 2, we can see 
that they are in better agreement (see Table ^). 

Table 4. Special episodes of superbolides: the new pres- 
sure of fragmentation calculated according to Eq. @. 
See Table 2 for other details. Pressures are expressed in 
MPa. 



Name 


Min P 


Max P 


S 


Pribram 


37 


110 


50 


Lost City 


6 


18 


50 


Sumava 


0.6 


1.7 


1 


Innisfree 


7 


22 


10 


Space based obs. 


8 


24 


50 


Space based obs. 


6 


18 


50 


Benesov 


2 


6 


10 


Peekskill 


4 


12 


30 


Marshall Isl. 


60 


180 


200 



Therefore, the new condition for fragmentation under 
unsteady regime is: 



V 



s 



k(1 + a)p s 



exp( 



h + H 
H 



(15) 



In addition, we have to consider that the distortion of 
the shock wave leads to the partial removal of the gas cap 
around the cosmic body, so that the ablation increases 
strongly. This enhances impact ionization, so that a — ► 1 
also during unsteady conditions, where thermal ionization 
is negligible. 

9. Examples 

Let us to consider two episodes in order to show how this 
"embryo" of theory works. Firstly, we can consider the 
Lugo bolide of 19 January 1993 (Cevolani et al 



1993 



Foschini 1998). It was a very bright bolide, which reached 
a peak magnitude of about —23 and released an esti- 
mated energy of about 14 kton, when exploded at about 
30 km over the city of Lugo, in northern Italy. In pre- 
vious analyses, it was considered that the fragmentation 
occurred when the dynamical pressure reached the me- 
chanical strength. 

If we now apply the condition given from Eq. (|i~5|), 
taking into account a — > 1, p s i = 1.293 kg/m 3 , and that 
the scale height H = 6.8 km at 30 km height, we obtain 
the values in Table |5|. 

We have three reasonable solutions, which are under- 
lined. If we consider that the final airburst occurred at 
about 30 km height, typical for type I bodies (see the dis- 
cussion at the end of Sect. 7), we can consider 26.9 km/s as 
the most probable speed. This value is in agreeme nt wit h 
first estimation from eyewitnesses (Cevolani et al. 1993] ). 
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Table 5. Values of the speed of the Lugo bolide at the 
moment of fragmentation [km/s]. 
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S [MPa] 


K = 2 


K = 6 


1 


6.6 


3.8 


10 


20.8 


12.0 


50 


46.5 


26.9 


200 


93.1 


53.7 



For the Tunguska event, a more detailed analysis is 
in preparation in collaboration with the members of the 
Tunguska99 Scientific Expedition. Here we want to under- 
line only one thing: from Fig. [I] we can see that, in subartic 
summer, the temperature does not change in the height in- 
terval crossing the troposphere. Therefore, the Mach num- 
ber does not change and we can apply the Eq. (|T2|), for 
steady state conditions, as sho wn in the Paper I. The er- 
ror, noted by Bronshten ( [2000 ), introduced in considering 
the explosion height, instead of the fragmentation height 
is negligible: the new value is 16 km/s to be compared 
with the old one of 16.5 km/s. 



10. Conclusion 

In this paper, we have done a further step toward the 
construction of a theory for the fragmentation of a small 
asteroid during the atmospheric entry. In the Paper I, we 
showed a specific part of the theory applied to a particular 
episode, the Tunguska event of 30 June 1908. 

Here we showed more details, both explaining better 
some assumptions in the theory of steady state motion 
outlined in the Paper I, and extending the theory to un- 
steady motion. We have taken into account the effect of 
turbulence and its interaction with shock waves. Some 
examples are discussed and we have found a reasonable 
agreement with available experimental data. 

We have found two conditions for the fragmentation, 
according to steadiness of the motion. 

1. Steady state motion: in this case, the compressibility 
suppresses the turbulence and, therefore, we can use 
Eq. @. 

2. Unsteady motion: we have a strong interaction be- 
tween the shock wave and the turbulence, which give 
rise to sudden pressure outburst. We have to use 
Eq. ©. 

On the other hand, it is necessary to remember that we 
are speculating, because of the scarce experimental data 
and the large uncertainties affecting records. Therefore, 
these researches must be taken cum grano salis. 
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